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ABSTRACT 

We examine the 3-10 keV EPIC spectra of Mrk 205 and Mrk 509 to investigate 
their Fe K features. The most significant feature in the spectra of both objects 
is an emission line at 6.4 keV. The spectra can be adequately modelled with a 
power law and a relatively narrow (a < 0.2 keV) Fe Ka emission line. Better fits 
are obtained when an additional Gaussian emission line, relativistic accretion-disk 
line, or Compton reflection from cold material, is added to the spectral model. 
We obtain similar goodness of fit for any of these three models, but the model 
including Compton reflection from cold material offers the simplest, physically 
self-consistent solution, because it only requires one reprocessing region. Thus the 
Fe K spectral features in Mrk 205 and Mrk 509 do not present strong evidence for 
reprocessing in the inner, relativistic parts of accretion disks. 
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1 INTRODUCTION 

X-ray observations probe the central regions of AGN. In 
the standard paradigm, this corresponds to the inner parts 
of an accretion disk around a supermassive black hole. 
Above the disk, a hot corona Compton upscatters optical- 
EUV photons to X-ray energies; some of this X-ray radi- 
ation is reprocessed in the surrounding material includ- 
ing the disk, giving rise to prominent Fe Ka lines. The 
broad, distorted velocity profile of Fe Ka emission, sug- 
gesting an accretion disk around a supermasive black hole, 
was first observed with ASCA in MCG -6-30-16 (Tanaka 
et al. 1995). This profile, with a sharp blue wing and a 
broad red tail, is a remarkable probe of the strong gravity 
regime. 

Studies of other AGN with ASCA suggested that 
broad, low-ionisation Fe Ka emission is common in AGN 
(Nandra et al. 1997, Reynolds et al. 1997), but it was not 
until the launch of XMM-Newton that the diversity of Fe 
line profiles could really be investigated. With the large 
increase in collecting area afforded by XMM-Newton, it 
was soon noticed that some luminous AGN showed differ- 
ent Fe line profiles to that of MCG -6-30-15 (e.g. Reeves 
et al. 2001a, Blustin et al. 2002). One particularly inter- 
esting example, the luminous Seyfert 1 galaxy Mrk 205, 
appeared to have a narrow, neutral Fe Ka line at 6.4 keV, 
accompanied by a broad line from He-like Fe (Reeves et 
al. 2001b). Reeves et al. argued that the ionised Fe line 
originates in the inner parts of an accretion disk, while 



the neutral Fe Ka line originates in a molecular torus, hy- 
pothesised to lie outside the broad line regions in AGN 
unification schemes (Antonucci 1993). A later observation 
of another Seyfert galaxy, Mrk 509, showed a very similar 
pair of narrow-neutral and broad-ionised Fe Ka emission 
lines (Pounds et al. 2001), demonstrating that this config- 
uration of Fe Ka profiles is not an isolated phenomenon 
in Mrk 205. 

In this paper we revisit the Fe K features in Mrk 205 
and Mrk 509 seen by XMM-Newton. By coadding the spec- 
tra from the three EPIC instruments we are able to max- 
imise the signal to noise per bin while properly sampling 
the EPIC spectral resolution around Fe K. The paper is 
laid out as follows. In Section we describe the obser- 
vations and data reduction, and we describe the spectral 
fitting in Section [3] The results are discussed in Section 
2] and we present our conclusions in Section |S] The Ap- 
pendix contains a description of the method employed to 
coadd spectra from the different EPIC instruments. 



2 OBSERVATIONS AND DATA REDUCTION 

Mrk 205 was observed with XMM-Newton on the 7th 
May 2000, and these data were presented by Reeves et 
al. (2001b). Several exposures were taken in both MOS 
and PN cameras, in full frame and large-window modes. 
Spectra of the source were extracted from circular regions 
of radius ~ 50" and background spectra were obtained 
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Table 1. XMM-Newton observations of Mrk 205 and Mrk 509 



Object 

Mrk 205 
Mrk 509 
Mrk 509 



Date 

7 May 2000 
25 November 2000 
20 April 2001 



Exposure 
(ks) 
49.0 
23.4 
30.6 



Count rate 
(count s _1 ) 

4.9 

26.7 

38.3 



from nearby source- free regions. All valid event patterns 
(singles, doubles, and triples) were selected in MOS, and 
only single events in PN. The spectra were combined using 
the procedure outlined in Appendix A. 

Mrk 509 has been observed twice with XMM-Newton. 
The first observation took place on 25th November 2000 
and the data were presented by Pounds et al. (2001). The 
second observation was performed on the 20th April 2001. 
In both observations, MOS and PN cameras were operated 
in small window mode. Source spectra were taken from a 
circular region of 40" — 50" radius, and background spec- 
tra were obtained from nearby regions free from bright 
sources. Single events were selected in MOS, and single 
and double events were selected in PN. The spectra were 
combined using the procedure outlined in the Appendix to 
produce one spectrum for each observation and one spec- 
trum which is a combination of the two. 



3 RESULTS 

The spectral fitting was performed with XSPEC. Only the 
rest- frame 3-10 keV energy range was used in the spectral 
fitting because we are primarily interested in the Fe K 
features. The broad emission lines reported by Reeves et 
al. (2001b) and Pounds et al. (2001) are only significant 
between 5 and 8 keV (see Fig. 4 of Pounds et al. 2001), 
so the 3-10 keV energy range allows a good measurement 
of the continuum on either side of Fe K, while excluding 
the noisier and less-well calibrated data at higher energies 
and the complex spectrum found at lower energies. We 
included the small effect of absorption from the Galaxy as 
a component in all our spectral modelling (Nh = 2.9 x 10 20 
cm" 2 towards Mrk 205 and Nh = 4.1 x 10 20 cm" 2 towards 
Mrk 509, Dickey & Lockman 1990). The results of the 
spectral fits are given in Table |5| 



3.1 Mrk 205 

We began by fitting a power law model. The counts spec- 
trum is shown in Fig. along with the power law model 
convolved with the instrument response. Like Reeves et 
al. (2001b) we find this model is a poor fit, and there are 
significant residuals around 6.4 keV. We therefore added 
a gaussian emission line, and obtained an acceptable fit, 
with a resolved line consistent with 6.4 keV Neutral Fe Ka. 
Although the fit was acceptable, residuals remained at ~ 7 
keV, and so we tried further fits including one of three ad- 
ditional model components that might plausibly account 
for these residuals and improve the fit: a gaussian emis- 
sion line at E > 6.5 keV, a relativistically-broadened ac- 
cretion disk line (Fabian et al. 1989; the "diskline" model 
in XSPEC), and Compton reflection from cold material 
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Figure 1. 3-10 keV EPIC spectrum of Mrk 205 and power law 
model (top panel) and data/model ratio for the five different 
model fits (lower panels). 



("pexrav" in XSPEC; Magdziarz and Zdziarski 1995). In 
the diskline model the energy of the line was constrained 
to lie between 6.7 and 6.9 keV, corresponding to He-like or 
H-like Fe (as proposed by Reeves et al. 2001b). In the re- 
flection model the inclination was fixed at 45°, and we as- 
sumed Solar elemental abundances (Anders and Grevesse 
1989). The addition of any of these three components re- 
sulted in a similar goodness of fit. The energies of the sec- 
ond Gaussian line and diskline components are found to 
be consistent with Fe XXVI, and in the case of the diskline 
the best fit was found for a line produced many Re from 
the black hole giving it a narrow profile, similar to that of 
the Gaussian line. We have also tried more complex mod- 
els, including two emission lines as well as reflection: the 
lowest value of x 2 / u ~ 138/135 was obtained when the 
second emission line is a Gaussian. However, according to 
the F-test, this is only a la improvement over the model 
including a single emission line and Compton reflection. 
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Figure 2. 3-10 keV EPIC spectrum of Mrk 509 from the first 
observation and power law model (top panel) and data/model 
ratio for the five different model fits (lower panels). 



Figure 3. 3-10 keV EPIC spectrum of Mrk 509 from the second 
observation and power law model (top panel) and data/model 
ratio for the five different model fits (lower panels). 



3.2 Mrk 509 first observation 

As for Mrk 205, we began with a power law model and ob- 
tained a completely unacceptable fit with large residuals 
around 6.4 keV (see Fig- El • Addition of a Gaussian line at 
~ 6.4 keV resulted in a much better, but still poor fit (re- 
jected with > 99.9% confidence). Adding a further Gaus- 
sian, a diskline or a reflection component improved the 
goodness of fit, which is marginally better when the third 
model component is a diskline rather than a reflection 
component or a second Gaussian line. However, whichever 
additional model component is included the model is still 
unacceptable at > 99% confidence. We have also tested 
a more complex model, including a second emission line 
as well as reflection, and obtain a best x 2 l v = 182/138 
for a diskline; this is only marginally better (2a accord- 
ing to the F-test) than the fit with reflection but without 
the second emission line. Strong residuals at 8.3 and 9.8 
keV together contribute ~ 30 to the \ 2 > an d are there- 
fore largely responsible for the poor model fit. We have 
no credible explanation for these residuals other than as 
statistical fluctuations. 



3.3 Mrk 509 second observation 

A power law model provided a poor fit to the spectrum 
(see Fig.[3Jl, and once again we found that the addition of 
a ~ 6.4 keV emission line resulted in a significantly im- 
proved (and quite acceptable) X 2 l v ~ 124/144. The addi- 
tion of an extra Gaussian, diskline or reflection component 
each improved the fit further, resulting in very a good fit. 
However, the addition of a second emission line to the re- 
flection model results in no further improvement in X 2 lv. 

3.4 Mrk 509 both observations combined 

The overall spectral shape of Mrk 509 is extremely sim- 
ilar in the two observations, and so to improve signal to 
noise we have coadded the data from both observations. 
A model containing only a power law is convincingly re- 
jected, but a power law and a single Gaussian provide an 
acceptable fit to the data. Adding an additional Gaussian, 
diskline or reflection component makes a significant im- 
provement to the x 2 l v ( a l > 99.9% significance according 
to the F-test), resulting in a good fit to the data. Finally, 
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Figure 4. 3-10 keV EPIC spectrum of Mrk 509 (both ob- 
servations coadded) and power law model (top panel), and 
data/model ratio for the five different model fits (lower pan- 
els). 



we have tried more complex models combining reflection 
with a second Gaussian or accretion-disc emission line. In 
this case the \ 2 l v is poorer than for the model including 
reflection without a second emission line. 



4 DISCUSSION 

The most significant spectral feature in all the 3-10 keV 
EPIC spectra of Mrk 205 and Mrk 509 is low-ionization 
Fc Ka emission at 6.4 keV. In both sources, improved fits 
are obtained when an additional emission component is 
included, peaking at slightly higher energy; Gaussian or 
relativistic emission lines or cold reflection are all plau- 
sible forms for this additional spectral component. For 
Mrk 509, the spectral modelling requires that the 6.4 keV 
line is broad, FWHM > 5000 km s _1 unless this additional 
component is included. More complex models, combining 
a second emission line with cold reflection, do not produce 
significantly better fits for any of the spectra than models 
with cold reflection but without a second emission line. 



The 6.4 keV emission is a signature of reprocessing by 
cold material, but where does this reprocessing take place? 
Both galaxies show negligible intrinsic absorption in soft 
X-rays, and hence the material responsible for the 6.4 keV 
emission must lie outside the line of sight to the contin- 
uum source. Possible locations include the accretion disk, 
the molecular torus favoured by AGN unification schemes 
(Antonucci 1993), and the (optical) broad line clouds. The 
molecular torus and the accretion disk represent Compton 
thick targets, and reprocessing at these locations will re- 
sult in a Compton reflection component with an edge at 
7.1 keV as well as Fe Ka line emission (Matt, Perola & Piro 
1991, George and Fabian 1991). However, the broad line 
clouds are expected to be Compton thin (Shields, Ferland 
& Peterson 1995) and therefore the Fe Ka line emission 
will not be accompanied by significant Compton scattered 
continuum. 

Simulations by Leahy and Creighton (1993, see also 
Yaqoob et al. 2001) show that if the broad-line clouds 
have column densities of 10 23 cm" 2 , they would need to 
cover 50% of the sky, as seen by the continuum source, 
to produce the Fe Ka line of 50 eV equivalent width in 
Mrk 509. To produce the 100 eV equivalent width line seen 
in Mrk 205, the broad line clouds would need to surround 
~ 100% of the central source. These covering fractions are 
much higher than the typical broad line region covering 
fractions deduced from the ultraviolet (10% - 25%, David- 
son & Netzer 1979, Goad & Koratkar 1998), suggesting 
that the broad line regions are probably not responsible 
for the majority of the Fe Ka photons. Furthermore in 
Mrk 509, even if the broad line region does produce the 
Fe Ka line, something else must contribute an additional 
broad spectral component at slightly higher energy, or else 
the velocity width of the Fe Ka line is inconsistent with 
the width of the optical lines, which have FWHM of only 
2270 km s _1 (Wandel, Peterson & Malkan 1999). 

On the other hand, if the Fe Ka line originates in 
Compton thick material, the equivalent width of the line 
suggests that this material intercepts ~ 30 - 60 per cent 
of the emitted radiation in Mrk 205, and about half as 
much in Mrk 509 (George & Fabian 1991). Therefore in 
both AGN the strength of the Fe Ka line alone suggests 
that a significant Compton reflection component should 
be present. When reflection is included in the fit we find 
that > 55% of the radiation is intercepted by the reflector 
in Mrk 205, and > 40% in Mrk 509, slightly higher than 
inferred from the Fe Ka line but not greatly so. Thus a 
Compton thick reprocessor that intercepts a significant 
fraction of the primary X-rays can account for all the Fe 
K features in Mrk 205 or Mrk 509. This reprocessor could 
be a distant molecular torus, or the accretion disk itself. 
However, the velocity width of the Fe Ka line implies that 
even if the reprocessor is the accretion disk, little of the 
reprocessed emission comes from the inner, relativistic, 
parts of the disk. 

So what evidence do we have for emission from highly 
ionized Fe, potentially in an accretion disk? For Mrk 205, 
the fits with a second emission line are as good as those in- 
cluding a reflection component. However the second emis- 
sion line is narrow when it is modeled as a Gaussian, and 
similarly the best fit accretion disk line is one in which 
the emission is dominated from material in the outer disk, 
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resulting in a narrow line. Hence even if the spectrum of 
Mrk 205 does include Fe XXV or Fe XXVI, there is no 
evidence for relativistic broadening. 

For Mrk 509, the best-fit second emission line is 
broad, whether it is modeled as a Gaussian or as a disk 
line, and the steep emissivity index implies that most of 
the emission would come from the inner part of the disk. 
The disk would have an inclination of between 20 and 40 
degrees, in agreement with the findings of Pounds et al. 
(2001). However, the \ 2 l v f° r the fits including an accre- 
tion disk line are only slightly better than the fits including 
Compton reflection for the first observation of Mrk 509, 
and are slightly poorer for the second observation; the 
combined spectrum is as well fit by either model. Hence 
although the features in the spectrum of Mrk 509 can be 
fit with a model including both a distant, cold, Compton- 
thin reprocessor and relativistically broadened emission 
from highly ionised Fe in an accretion disk, they can be fit 
equally well by reprocessing and Compton reflection from 
distant, cold material. 

Thus in both Mrk 205 and Mrk 509, we find that the 
Fe K features can be explained by a single phase Comp- 
ton thick cold reflector. While the presence of reflection 
from the highly ionized, high velocity, inner parts of an 
accretion disk is not ruled out by these data, it is not 
unambiguously detected. 



5 CONCLUSIONS 

We have analysed the 3-10 keV XMM-Newton EPIC spec- 
tra of Mrk 205 and Mrk 509 to investigate the Fe K fea- 
tures in these objects. Acceptable fits can be obtained for 
models containing nothing more than a power law and an 
emission line at 6.4 keV, consistent with cold Fe Ka. How- 
ever, better fits are obtained when an additional spectral 
component is included in the model, either Compton re- 
flection from cold material or an emission line from ionised 
Fe; the goodness of fit is similar whichever component is 
added. In Mrk 205, there is no evidence for relativistic 
broadening of any emission line, but in Mrk 509 the best 
fit parameters for an ionised Fe emission line suggest that 
it might originate in the inner regions of an accretion disk. 
However, illumination of distant, cold material provides 
a simpler, self consistent explanation of the spectral fea- 
tures than models including reflection from highly ionized, 
relativistic material. Therefore, contrary to Pounds et al. 
(2001) and Reeves et al. (2001b), we do not find strong 
evidence in either object for reprocessing in the highly 
ionised inner parts of an accretion disk. 
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APPENDIX A: METHOD TO COMBINE 
SPECTRA FROM EPIC MOS AND EPIC PN 
CAMERAS 

The following method was designed with pulse height 
spectra and response matrices in standard 'OGIP' 1 for- 
mat in mind. We use the term 'response matrix' to refer to 
the product of the effective area and the energy^channcl 
redistribution matrix. 

We wish to combine a number of individual 
source+background pulse-height spectra to a single pulse- 
height spectrum. We label each individual spectrum with 
the index s — 1, 2, ...N spec , and we use the index i = 
1,2, ...Nchanis) to denote the channels of each spectrum. 
Each channel is assigned a 'nominal' energy range of 
ENOM min (i) < E < ENOM max (i). The number of pho- 
tons in each channel i of spectrum s is C(s,i). In gen- 
eral, each source+background spectrum C(s, i) has a cor- 
responding background spectrum B(s,i), with a scaling 
factor F(s) relating the geometric area and/or exposure 
times of the two spectra. Each of the original spectra has 
a corresponding response matrix, whose elements contain 



1 http://heasarc.gsfc.nasa.gov/docs/heasarc/ofwg/ 
docs/spectra/ogip.92_007/ogip_92_007.html 
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the effective area for a given channel and a given en- 
ergy range. The element of the response matrix for spec- 
trum s corresponding to a particular channel i and a par- 
ticular energy range E mi „(j) < E < E max (j) (where 
j = 1,2, ■■■N ranae (s)) is denoted R(s,i,j). Throughout we 
use capitalised indices when refering to the summed spec- 
trum and response matrix; the index s is omitted when 
referring to a combined spectrum or response matrix. 

We define the fractional overlap of channel i of spec- 
trum s with channel / of the combined spectrum to be 

if ENOM mm (s,i) < ENOM m ax (I) 
and ENOM max (s,i) > ENOM mm (I) 

f(s,i,I) = 

{min[ENOM max (s,i),ENOM max (I)]- 
max[ENOM min {s,i),ENOM min (I)}} / 
{ENOM max (s,i) - ENOM min (s,i)} 
otherwise 

/(«,*,/) =0 (Al) 

The combined spectrum is constructed by summing all the 
counts from all the spectra in the nominal energy range of 
each channel. So, 

"spec N chan (s) 

C W= E E /(M,/)C(M) (A2) 

s=l i=l 

If the nominal energy ranges of the channels of the dif- 
ferent spectra do not coincide with those of the combined 
spectrum, some randomization of photons will be required 
to ensure that the channels of the output spectrum contain 
integer numbers of counts. 

We combine the individual background spectra into 
a single background spectrum with a single scaling factor 
F. The following scheme produces a background photon 
spectrum with the signal to noise ratio propogated from 
the individual background spectra. 

N spec N chan (s) 

B ^ = T E F ^ E B ^ (A 3 ) 



if E m in(s,j) < E max (J) 

and E max (s,j) > E min (J) 

g(s,j, J) = 

min[E max (s,j), E max (J)] - max [E mln (s,j), E min (J)] 
otherwise 

g(a,j,J)=0 (A5) 
The combined response matrix is then constructed accord- 
ing to: 
R(I,J) = 

" v -"""'" E":r 9e{s) 9(s,j,j)R(s,i,j) 



E E /(^ J )- 



(A6) 



where 



B(s,i) 



Ef= r F{s 



(A4) 



i£tr n(3) B(s,i) 

The response matrices are easily combined without 
any complicated weighting, provided that the pulse height 
spectra are all realisations of a single spectrum (e.g. if 
they are from observations of the same source at the 
same time but by different instruments). If the individ- 
ual pulse height spectra are realisations of intrinsically 
different spectra, then the response matrix combination 
described here will not generally be appropriate. However 
if the spectra differ only in intensity then a simple scaling 
factor can be used to weight the contributions of the dif- 
ferent response matrices to the final spectrum. For each 
energy range of the combined response matrix, the ele- 
ments corresponding to each channel are combined in the 
same way as the original spectra. This is more complicated 
if the energy ranges differ between the response matrices. 
In this case we use a weighted average response for a given 
energy range. We define the fractional overlap of energy 
range j of spectrum s with the energy range J of the out- 
put response matrix as follows: 
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